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I. INTRODUCTION 

The low dimensional electron and spin systems are usu- 
ally regarded as important playgrounds for the studies 
of correlated quantum matters due to the strong quan- 
tum fluctuation effects in one and two spatial dimensions. 
Many exotic non-Fermi liquid states such as the Lut- 
tinger liquid and spin liquid states, to name just a few, 
are found in the single or coupled chain systems. How- 
ever, nature never stops surprising us. The recent ad- 
vances in the technique of loading ultracold gases into op- 
tical lattices open a new era in the research of condensed 
matter systems. Among the recent achievements, the ex- 
perimental realization of the Boson-Hubbard model sig- 
natures an important step in this direction, not only 
because of the tunability of the controlling parameters 
in the corresponding experiments, but also because it 
facilitates the first observation of the Mott insulating 
state of bosons and the associated quantum phase tran- 
sitioni While the atomic interactions in the ultracold 
gases can be treated as contact ones for most cases, a siz- 
able longer range interaction is now within experimental 
reach by using the dipolar interaction among atomS)^^ 
which could provide further opportunities of control- 
ling/designing new experiments. Following these lines 
of developments, a natural question to ask is whether 
or not the longer range interactions can trigger a stable 
new quantum phase of matter which contains non-trivial 
internal structure. 

Among many research works devoted to the under- 
standing of the effects of long-range dipolar interactions, 
we mention the recent work by Emanuele G. Dalla et 
al^ who studied the one-dimensional boson insulators 
within the context of an extended boson Hubbard model 
(EBHM) by employing the density matrix rcnormaliza- 
tion group (DMRG) method. By tuning the ratio of the 
on-site interaction over the hopping amplitude U/t and 
the ratio of the longer range interaction over the hopping 
amplitude V/t, the mean field analysis shows that three 
different conventional phases can be reached. These in- 
clude the Mott insulator at large U, the density wave 
state for large V and a superfluid state for large t. The 



surprising thing is that a new intermediate insulating 
state, the Haldane insulator, which separates itself from 
the other two insulating states by second order quantum 
phase transitions was found in Ref. [U. Moreover, it was 
shown that such a state possesses a non-vanishing non- 
local string order, similar to the Haldane phase of the 
quantum spin-one chain. Such a state is definitely be- 
yond the reach of the traditional one-dimensional hydro- 
dynamic effective theory for the one-dimensional boson 
superfluid-to-Mott transition^ Recently, the present au- 
thor and his collaborators have developed an phenomeno- 
logical two- component hydrodynamical effective theory^ 
which successfully captures the main features of all the 
phases, including the Haldane insulator, found in the re- 
cent DMRG study of the one-dimensional EBHM at in- 
teger fillings. This effective theory also clarifies the na- 
ture of the quantum phase transitions between different 
phases. 



Knowing the above results, it is desirable to see if sim- 
ilar exotic phases or insulating states can be found in 
other one-dimensional or quasi-one-dimensional systems. 
As a first attempt, we consider in this paper a dipolar 
boson system of two weakly coupled chains within the 
framework of the EBHM. Interestingly, we found that 
the competition of the inter-chain hopping and the inter- 
chain interaction does lead to two different types of Mott 
insulating states, with one of them possessing a nontriv- 
ial string order. In addition to that, we also found that 
the inter-chain attraction can give rise to an interesting 
paired superfluid state where the inter-chain bound bo- 
son pairs show an algebraic long range superfluid order 
while the single-boson superfluid correlations decay expo- 
nentially The rest of the paper is organized as follows: 
In section II, we introduce our model and its effective 
theory. In section III and IV, we analyze the phases of 
the model and discuss the issue of the string order in 
the Mott insulating state. The final section is dedicated 
to our conclusions, and the resulting phase diagram is 
summarized in figure 2. 
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FIG. 1: (color online) The two-chain lattice model described 
by Hamiltonian The arrows denote the possible orienta- 
tions of the dipole moments. The left figure (a) corresponds 
to the IDW phase and the right figure (b) represents the BDW 
phase considered in the paper. (See section III.) 



II. THE MODEL AND ITS EFFECTIVE 
THEORY 

The lattice bosons studied in the present paper is de- 
scribed by the following extended Bose-Hubbard model 



H = -t E (^+H.c.) + ^E^ 



V ^2 8n ai 8i 



<*3 



+ t±_ E( fot ii 6 2i + h.c.) + V± E SnuSnai , (1) 

i i 

where a = 1,2 is the chain index, t is the nearest- 
neighbor hopping along the chain, t± is the inter-chain 
hopping, U is the on-site repulsion, and V, V± arc 
the nearest-neighbor interactions along and between the 
chains, respectively. The operator b* ai create a boson at 
site i on chain a, n a i = b^bd is the number operator, 
and drio-i = n a i — n measures the deviation of the particle 
number from a mean filling n. In the present paper, we 
will focus our attention to the cases of half-odd-integer 
fillings, i.e. n = N + 1/2 where N is a non-negative 
integer. 

We stress that the above model can be realized by po- 
lar molecules,— £ or atoms with a larger dipolar magnetic 
moment such as 53 Cr4 Moreover, it is possible to adjust 
the mutual orientations between the dipoles by using ex- 
ternal electric (magnetic) fields, so that we may control 
the signs of the nearest-neighbor interactions V and V±. 
Therefore, both the cases of the attractive and the repul- 
sive nearest-neighbor couplings will be considered below. 
(See Fig. [TJ) With this understanding in mind, we now 
discuss the low energy effective theory of the model de- 
fined in Eq. Q. 

For simplicity, we first consider the case where the on- 
site repulsion U is the largest energy scale, i.e. U » 
t,V,t±,V±. Under this condition and near half-odd- 
integer fillings, we can truncate the boson Hilbcrt space 
so that only states with local occupations n a i = N+l and 
N are allowed, and the boson operators can be mapped to 



spin-1/2 operators, b\ i — > 5+, 5n ai — > S z oi . In this low- 
energy subspace, the extended boson Hubbard model is 
mapped to coupled xxz spin-1/2 models, 

H = H + H ± 

Fo = -t E ( S ?i 3 ij + h - c -) + V E S 'i S U + C 1 - 2 ) 

(»>i) <«>i> 

ff_L = t ± E(^i^ + h-c) + V ± J2 S u S 2i ■ (2) 

i i 

In the following, we will assume that the inter-chain 
coupling tj_ , Vj_ are much smaller than the intra-chain 
couplings t and V. Then the spin-1/2 chain can be 
bosonized using the standard method^ and the intra- 
chain Hamiltonian can be written in the following form, 



x ^{d x ^f + K{d x 6„Y 



I dx COS V ^4>cr ■ (3) 

At the perturbative level, t 3> V, we have v = 
v \Jl + ^ = ^,v = too, and g = (a is the 

lattice spacing.). However, the validity of this effective 
action goes beyond its perturbative derivation, and the 
relations between the Luttinger parameter K and v and 
the spin-chain couplings t, V can be established exactly, 
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with A = Yt m the above equation. 

The inclusion of the inter-chain coupling terms do not 
present any new difficulty. Upon including them, the 
resulting Hamiltonian becomes 



*= E 
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In the above equation, <fi s , 9 S and <p a , 9 a are the symmetric 
and anti-symmetric combinations of the boson fields for 
the spin-1/2 operators, respectively. In terms of them, 
the original lattice boson operators can be expressed as, 



h/2,i 
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and the bond and inter-chain density fluctuations are 
Sn 
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Since the intra-chain nonlinear term 
g cos \/&K(j) s cos ^/%TT(f)a is always less relevant than 
the inter-chain coupling terms in the region that we are 
interested in below, it can be safely neglected. Therefore, 
in the strong on-site repulsion limit, the two-chain ex- 
tended boson Hubbard model at half-odd-intcger fillings 
can be described by the following effective hydrodynamic 
theory, 
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, = , and v,K here are pa- 

ramctcrs defined in Eq. (|4|). At this point, it is interesting 
to notice that exactly an action of the same form occurs 
in our recent study of the one-dimensional boson Mott 
transition near integer fillings^ albeit in a completely 
different physical context. 

Before we analyze the consequences of the Hamiltonian 
in Eq. ([5]) , we emphasize that this effective action is quite 
general and is not restricted to hard-core bosons. For the 
finite U soft-core bosons, one may introduce the boson 
field and density operators,— 

L. = f po + * e iM*pa*+V*<i>~) e W*9* 



Of) 



(9) 



where a = 1,2 are the chain indices, and po f» (N+i)/ao 
is the average boson density. After some simple manip- 
ulations, it is not hard to see that we will still get the 
same effective action Consequently, the results we 
get in this paper is robust as long as the inter-chain cou- 
plings are weak. The only problem is that in soft-core 
boson case, the relations between the Luttinger param- 
eters K sa and the microscopic couplings are unknown 
and Eq. can only be treated in a phenomcnological 



III. THE PHASES AND PHASE TRANSITIONS 
OF THE EXTENDED BOSON HUBBARD 
MODEL 

After establishing our effective theory for the two-chain 
EBHM, we now turn to study the possible phases and the 
associated quantum phase transitions of this system. The 
first thing to notice is that the charge U(l)xU(l) symme- 
try of the decoupled chains is broken down to a U(l)xZ2 



by the inter-chain hopping term t±(b\ i b2i + h.c), which 
can be viewed as the inter-chain Josephson coupling be- 
tween two one-dimensional boson superfluids. We also 
notice that the diagonal U(l) symmetry here forbids 
terms like cos (39 s . In addition to that, there is a lattice 
translation symmetry of the EBHM. As we shall see im- 
mediately, these discrete symmetries can be broken spon- 
taneously due to correlation effects. 

The only nonlinear term which can potentially open a 
gap in the symmetric mode is the gs term in Eq. l[8]). 
and its dimension is 2K S . Therefore, we expect that 
depending on whether K s > 1 or K s < I, the sym- 
metric hydrodynamic mode will either be gapless or ac- 
quire a gap. Furthermore, the transition correspond- 
ing to the opening of the spectral gap in the symmetric 
mode falls into the KT universality class. On the other 
hand, for the antisymmetric mode, there are two com- 
peting nonlinear terms gi and 172 in Eq. l[8]). which have 
scaling dimensions ■^k- and 2K a , respectively. Hence, 
depending on whether K a > 1/2 or K a < 1/2, either 
the gi cos \[2~KB a term or the 32 cos ^/8TT(f> a term will be 
dominant and opens a gap in the anti-symmetric mode. 
In other words, the anti-symmetric mode will always be 
gapped except at the critical line defined by K a = 1/2, 
where the critical theory falls into the Ising universal- 
ity class. The nature of this transition and the asso- 
ciated critical mode can be seen as follows^ Exactly 
at K a = 1/2 ji£ we may introduce a set of right- and 



left-moving Dirac fcrmions xr,l = 



i e ±V^<P' RiL with 



/'2iva 



T\/ T H « = ®a + T\j\Qa to rewrite 
In order to identify the critical mode 



VR,L ~ - V 2K a 

the Hamiltonian. 
and the physical low lying excitations, it is convenient to 
further decompose the Dirac fcrmions into their real and 
imaginary parts Xv = ^7 2 = (^ + *P^)> v — R,L, and express 
the Hamiltonian in terms of these Majorana fcrmions 



H ^(prOxPr- PlOxPl) +im p p R p L . (10) 



In the above equation, the Majorana fermion masses are 

™e = $f . m p = when Isal = \9z\, one of 

the Majorana fermion will become gapless and the re- 
maining one will still be massive. Hence, we see that the 
competition between the inter-chain tunnelling and the 
inter-chain interaction may result in an Ising transition. 
The above discussions suggest that, depending on K s 1 
and K a 1/2, there can be four different phases. In the 
following, we shall elaborate on the natures of these four 
phases. 

We start by considering the case where K s > 1 and 
K a > 1/2. In this case, the only relevant part of the 
interaction Hamiltonian is the gi-term, which arises from 
the inter-chain Josephson coupling. Using Eq. ([5]) and 
Eq. ([7]), one may easily see that the single particle and 
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density correlations are 



{8n ai 5n a j) 
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an exponentially decaying 
Vstaggerred part 



(11) 



It is also interesting to examine the following long dis- 
tance two-particle correlations, 



\i-j\l/K. 



(bl&iblfoj) ~ const. 



(12) 
(13) 



From the above equations, we see that the single-particle 
correlation is dominant over the two-particle pair corre- 
lation, as is expected for a usual superfluid. However, 
Eq. (fl~3|) also shows that the system exhibits a true long- 
range order for the pair condensate (b' u b2i) 7^ 0. The ex- 
istence of a non-zero condensate in this one-dimensional 
system indicates that some kind of resonant bond-pairs 
are formed for bosons on two different chains, and it is 
a direct manifestation of the phase-locking between the 
two chains due to the inter-chain Josephson coupling. 
We also notice that although the density correlations 
along the chain decay algebraically, the density-difference 
5n_i =: V u bii — b\^>ii '■ nas an exponentially decaying 
correlation function. This is consistent with the above 
resonant bonding-boson pair picture. Since there is no 
translational symmetry breaking in this phase, and mo- 
tivated by Eq. lfTTj) and Eq. (fT3|) . this phase will be coined 
as the single-particle resonant superfluid (RSF) later in 
this paper. 

We next consider the second gapless phase specified by 
the Luttinger parameters K s > 1 and K a < 1/2. In this 
case, the relevant part of the interacting Hamiltonian is 
(?2-term which originates from the inter-chain dipolar in- 
teraction between bosons. The single particle and density 
correlations in this phase are 



(b\ib\j) = {b\jb2j) ~ decays exponentially , 

1 f-lV^'l 
+ {const.) x ). J . IK ^ . (14) 
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The most interesting feature of this phase is that al- 
though the symmetric mode is gapless and possess a non- 
zero superfluid stiffness, the single-particle superfluid cor- 
relation decays exponentially. On the other hand, two- 
particle correlation shows an algebraic long range order, 



1 



|i-i|V*. ' 



decays exponentially 



(15) 



From the above results, we conclude that this transla- 
tionally invariant phase is a one-dimensional paired su- 
perfluid (PSF), where the inter-chain bond-boson pairs 
flow coherently along the chain. 



We now turn to discuss the remaining two gapped 
phases. For K s < 1 and K a > 1/2, both the inter- 
chain Josephson and the inter-chain dipolar interaction 
are relevant, and these interactions open gapes for both 
the symmetric and ant-symmetric modes. The relevant 
part of the effective Hamiltonian in the present case is, 



H = 



E y 1 



-^(d x ^ a ) 2 + K s (d x 6 a ) 2 



+ .9i J dx cos \phnQ a + 53 / dx cos "\Z87r0s • (16) 

In this phase, it is easy to see that among the correlation 
functions examined previously, the only one which does 
not decay exponentially is the resonant-pair condensate 
Eq. (|13|) . This is a translationally invariant insulating 
state with a non-zero resonant pair boson condensate. 
Therefore, it is a Mott insulator (MI). 

The last phase corresponds to K s < l,K a < 1/2, and 
the relevant part of the effective Hamiltonian in this case 



H 



T 

a— s.a 



dx 



K, 



■{d^f + K s {d x e a f 



+ 92 / dx cos v 87T0 a + <?3 / dx cosv87r</> s , (17) 



with —33 = 52 oc V± in the simplest EBHM. In this 
phase, all the single-particle and the pair correlations de- 
cay exponentially. On the other hand, for the bond- and 
inter-chain density fluctuations, 8n±i cx: b^bu ± b\^>2i '■, 
we have 



V x > 0, (Sn. 
Vx < 0, (Sn, 



(-l) x '/a Q 



(18) 



This result indicates that the discrete Z2 lattice transla- 
tion symmetry is broken spontaneously. Therefore, de- 
pending on the sign of the inter-chain dipolar interaction, 
the system exhibits either a bond-density wave pattern 
(BDW) or an inter-chain density wave pattern (IDW). 
(See Fig.l for a pictorial illustration.) Notice that if we 
change the relative orientations of the dipolcs between 
the two chains so that the inter-chain dipolar interaction 
changes its sign, we will tunc a quantum phase transition 
across these two different density wave phases and the as- 
sociated critical theory is a U(l)xU(l) Gaussian theory. 
These translationally non-invariant states can, in prin- 
ciple, be observed experimentally in the corresponding 
spatial image patterns of the time-of- flight experiments. 



IV. THE HIDDEN STRING ORDER 

Inside the Mott insulator phase, we can also adjust 
the relative orientations of the dipoles between the two 
chains to change the sign of Vx- This will introduce a 
further quantum phase transition due to the accidental 
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FIG. 2: (color online) The phase diagram of the two weakly 
coupled half-filled one-dimensional EBHM obtained in the 
hard-core boson limit. Notice that for the hard-core bosons, 
only the BDW phase can occur. However, we expect that this 
is a artifact of the hard-core condition, and for more general 
cases, depending on the sign of V±, both the IDW and the 
BDW phases can be stablized. See the discussions in the last 
section. Notice also that the phase transition across the red 
line is of the Ising universality class, while the phase transition 
across the blue line belongs to the KT universality class. The 
Mott insulator II below the green line contains a non-trivial 
string order, and it is separated from the Mott insulating I 
state by a U(l) Gaussian critical theory. 

gap closing in the symmetric mode alone, and the corre- 
sponding critical theory is a U(l) Gaussian theory. On 
both sides of the transition, the system is gapped and 
do not break any lattice symmetry. So, what is the dif- 
ference between these two Mott insulating states? Inter- 
estingly, from the structure of the above hydrodynamic 
theory and the results of recent studies on the ID boson 
Mott insulator, — £ we expect that one of the Mott insu- 
lating phase will posses a non-vanishing non-local string 
order, — while the other does not. This can also be seen 
by noticing that our effective Hamiltonian bears some 
similarities with that of isotropic spin ladders ; 12 ' 13 where 
nontrivial string orders were also found. More specifi- 
cally, drawing the analogy with the isotropic spin ladders, 
we may introduce the string operator 

d string (i-j) = lim (5n +l e"^=^ Sn ^5n +3 ). (19) 

\x-y\—too 

Following of the discussions of M. Nakamura^ we expect 
that the correct bosonized form of the string operator 
should be 

O string ~ lim ( sin V2n(j) s (x) sin V2n(j) s (y) ) , (20) 

\x— y\ — >oo \ / 

which is nonzero for V± < 0, due to the pining of the 
symmetric mode at <f> s = \p$- On the other hand, the 
string order vanishes for V±_ > where the field 4> s is 
pinned at (j) s = 0. Therefore, for the negative V± case, the 
MI phase we found here is the boson analogy of the AKLT 
phase of a spin ladder system, while for the positive V± 
case, the boson MI in the present EBHM is the boson 
analogy of the RVB phase of the spin ladder^ 



A more interesting question is perhaps about how to 
detect the differences between the two Mott insulating 
states. From the structure of the effective theory, the 
bulk spectrums of these two insulating states in the im- 
mediate neighborhoods of the transition line are almost 
identical. On the other hand, it can be shown that the 
MI state with a non-vanishing string order contains a 
zero energy edge stated This fact is most easily seen 
by examining the Sine-Gordon model of the symmetric 
mode near its Luther-Emery point^ (K s = 1/2). Near 
the Luther-Emergy point, the symmetric mode can be 
mapped to a massive Dirac fermion theory^ 

/>oo 

H= dx¥(x) [-ivo- 3 d x + cr' 2 m] + • • • , (21) 

Jo 

where ^(x) = (ip R (x),tpL(x)) T is a Dirac fermion field, 
m = —Vj_/2, and • • • denotes the possible four fermion 
interactions. These fermions should be identified as the 
low energy elementary (soliton) excitations of the sym- 
metric mode in the insulating states. For a system with a 
boundary at, say x = 0, we may introduce the boundary 
condition iPr(0) = - 0L(O)r^ and reformulate the above 
Dirac theory in terms of a single chiral component by 
introducing ipn(x) = iPl(—x) for x < 0. After doing so, 
the Dirac Hamiltonian reduces to 

/oo 
dx - ivipR(x)^ d x ipR(x) + imsgn(x)il> R (x)%l) R (-x) . 
-oo 

(22) 

It is not hard to see that the corresponding Schrodingcr 
equation 

- ivd x tpR(x) + imsgn(x)-ip R (-x) = eip R {x) , (23) 

contains a zero energy state localized at the edge, 

4> Re {x) ~ e~ m ^ , e e = 0. (24) 

Notice that this edge state exists only for m > 0, i.e. for 
V± < 0. Since all other bulk excitations have finite gapes, 
the existence of this zero-energy edge state constitutes a 
unique feature of the string-ordered MI state, and its ex- 
istence can be revealed in a zero frequency delta-function 
like peak in the optical absorption spectrum. Since in real 
experiments, the system is trapped in a finite domain, it 
is interesting to see if such a delta-function absorption 
peak can be observed experimentally. 

V. CONCLUSIONS AND DISCUSSIONS 

Although most of our results are obtained in the weak 
inter-chain coupling limit, we point out here that some 
of the important features of the weak coupling phase dia- 
grams can be realized even in the strong inter-chain cou- 
pling limit. For example, the effect of a large V± can be 
understood by considering the hard-core boson limit and 
keeping only the states with local occupations n a i = 0, 1. 
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Under this condition, we may introduce the bond opera- 
tors-^ 



*t| ) 



V2 



(|io> - |oi», 4|o) 
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4 |o> = ^=(|ii) + |oo)) , 410) = -W|io) + |oi)) , 
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V2' 



(25) 



nian becomes 



Heff = ~t± b ifi+> ~ tJ2( b +i b +,i+l + h - C -) 

i i 

U ^2 ^V- 

+ T nb+l : ' + J : nb+l :: Ub+: > : 



), 



(28) 



for each site, and the ket, say |10), denotes a state with a 
local bond-boson occupation \nu = 1, rin = 0), and these 
operators obey the constraint sjs, + ^2 a t ai t a i = 1. In 
terms of these bond operators, the inter-chain interaction 
can be written as 

H v± = vx $>L*x< + - + 4**i) ■ ( 26 ) 



When \V±\ S> \t\, \t±\, and Vj_ > 0, we may keep only 
Si,t z i, and the leading order effective low energy Haml- 
tonian becomes 



V 



v v — > 

H = -Vx + ^(nAiUAi+i + n B in Bl+1 ) 



V 



t 2 A 



^2(nAin B i+i + nBWAi+i) + 0( ^ ) , (27) 



where A\ = ^{s\ + tl),B} = ^(s\ - 4) create 
bond-pairs |10)i and |01)i, respectively, and riA.Bi = 
A\Ai, B\Bi are the number operators of the A and B 
bosons, respectively. We see from this strong coupling 
analysis that in the large Vj_ limit, the ground state con- 
figuration is completely determined by minimizing the 
inter-particle interactions between the A and B parti- 
cles. Therefore, for V > 0, we expect a Wigner crystal 
like state with the A and B particles are arranged in 
an alternating pattern • ■ ■ ABAB ■ ■ ■ , while for V < 0, 
phase separation will occur. In terms of the original 
two-chain boson language, this Wigner crystal like den- 
sity wave state is exactly the IDW phase discussed pre- 
viously. When Vx < 0, we get similar results, ex- 
cept now that the A and B operators are defined by 

A = ^(4 - tt\ilB\ = ^(4 + *'tji)' which crcatc 

states with local-bond occupations |ll)j and |00)j, re- 
spectively. Again, for V > the ground state is a density 
wave state and for V < we have phase separations. It 
is not hard to see that the Wigner crystal like state in 
this situation is just the BDW phase discussed earlier. 

Similarly, when the inter-chain hopping amplitude t± 
is the largest energy scale, we may introduce the inter- 
chain bonding and anti-bonding bosons fields b±^ = 
^|(&ij±&2i)- At low energy, the anti-bonding bosons can 
be projected out and the resulting low energy Hamilto- 



to the leading order, where Hb + i = b^ +i b + i is the density 
operator of the bonding bosons. This effective Hamilto- 
nian is nothing but that of the one-dimensional EBHM of 
the bonding-bosons at integer fillings. Transforming back 
to the original two-chain boson language, it is not hard 
to see that the supcrfluid phase of the bonding boson at 
large t, where the boson density distributed uniformly, 
automatically possesses a non-vanishing resonant-pair 
condensate (& li &2i + h.c). Hence this phase may be iden- 
tified as the RSF phase of the two-chain EBHM consid- 
ered in this paper. Similarly, its density wave phase at 
large V can be viewed as the BDW phase of the original 
two-chain system. According to Ref. [a, when V increases 
to the order of U, we expect the bonding boson will enter 
into a "Haldanc insulator" phase which has a non-trivial 
string order. Since in the truncated Hilbert space where 
the anti-bonding bosons are projected out, the density 
operator n,b + i coincides with the total bond boson num- 
ber operator n+j, Hence, the string operator defined in 
Eq. (fT9)) reduces to that of Ref. [H in the large t± limit. 
Although the transition between the Haldane insulating 
phase and the conventional Mott insulating phase of the 
bonding-bosons is achieved by tuning the intra-chain cou- 
pling V, while it is the inter-chain coupling Vx which 
plays a decisive role in the corresponding phase transi- 
tion of the weakly coupled chain case, the fact that the 
string operators coincide in certain limit suggests that 
these two phases are, in fact, adiabatically connected. 

We now briefly discuss some earlier results which are 
related to the present work. A phase analogous to our 
PSF state was noticed in the context of a coupled spin 
modcl^ and it was also shown to appear in a boson 
ladder model at incommensurate fillings^ The Mott- 
supcrfluid transition for a boson ladder model, coupled 
by the inter-chain josephson coupling alone, at the com- 
mensurate filling of one boson per site was studied in 
Ref. [H, while in our present case, the Mott insulating 
state occurs at half-odd-intcgcr fillings. We also notice 
that the hidden order in single chain electronic Hubbard 
models was studied using both the bosonization and the 
DMRG methods in Ref. [23|, and there arc some recent 
efforts which address the issue of string orders in elec- 
tronic band insulators and their relations with the Hal- 
dane chain^ These works, together with the results of 
Ref. and our work, suggest that the idea of classify- 
ing the boson or fcrmion Mott insulators in quasi-one- 
dimensional systems using the string order may be guite 
general. 
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To summarize, the above strong coupling analysis 
shows that the phases we obtained, when one of the 
inter-chain coupling is dominant, is consistent with the 
results we gained in the weak-coupling bosonization anal- 
ysis. Since one of the most important conclusion that 
follows from the bosonization study is that the interplay 
between the inter-chain hopping and inter-chain interac- 
tion results in two exotic phases — the paired supcrfluid 
state and the Mott insulating state with a string order, 
it is tempting to speculate that the major results of the 
bosonization analysis in this paper should be valid even 
when the inter-chain coupling strength is not so weak 
and can be extended to larger regions in the whole phase 
diagram. 

Using the analytic results obtained in this paper as 



a guide, we hope that it is helpful for future numerical 
works to determine the exact phase boundaries of the 
two-chain EBHM, and more importantly, the exotic in- 
sulating phase and the paired superfluid phase can be 
observed in the future experiments. 
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